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1. Introduction 



X. 

The generalized exponential-cosine-screened Coulomb (GECSC) potential or the generalized 
cosine Yukawa (GCY) potential: 

V (r) = — (—) exp(— 5r) cos(gdr), (1) 



r / 

where A is the strength coupling constant and 5 is the screening parameter, is known to 
describe adequately the effective interaction in many-body enviroment of a variety of fields 
such as atomic, nuclear, solid-state, plasma physics and quantum field theory [1,2]. It is 
also used in describing the potential between an ionized impurity and an electron in a metal 
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[3,4] or a semiconductor [5] and the electron-positron interaction in a positronium atom in 
a solid [6]. The potential in (1) with g — 1 is known as a cosine-screened Coulomb potential. 
The static screened Coulomb (SSC) potential (g — case) is well represented by Yukawa 
form: V(r) = —(aZe 2 )exp(—5r)/r which emerges as a special case of the ECSC potential 
in (1) with A = aZe 2 , where a = (137.037)" 1 is the fine-structure constant and Z is the 
atomic number, is often used for the description of the energy levels of light to heavy neutral 
atoms [7]. It is known that SSC potential yields reasonable results only for the innermost 
states when Z is large. However, for the outermost and middle atomic states, it gives rather 
poor results. Although the bound state energies for the SSC potential with Z — 1 have been 
studied [7]. 

The Schrodinger equation for such a potential does not admit exact solutions. So various 
approximate methods [8] both numerical and analytical have been developed Hence, the 
bound-state energies of the ECSC potential were first calculated for the Is state using 
numerical [3,8,9] and analytical [10,11] methods and for the s states by a variational method 
[12]. Additionally, the energy eigenvalues of the ECSC potential [13] have been recalculated 
for the Is state with the use of the ground-state logarithmic perturbation theory [14,15] and 
the Pade approximant method. The problem of determining the critical screening parameter 
S c for the s states was also studied [16]. 

It has also been shown that the problem of screened Coulomb potentials can be solved 
to a very high accuracy [17] by using the hypervirial relations [18,19,20] and the Pade ap- 
proximant method. The bound-state energies of the ECSC potential for all eigenstates were 
accurately determined within the framework of the hypervirial Pade scheme [21]. Further, 
the large-N expansion method of Mlodinow and Shatz [22] was also applied to obtain the 
energies of the ground and first excited s— states and the corresponding wave functions. 
Recently, we studied the bound-states of the ECSC potential for all states using the shifted 
large N— expansion technique [23]. 

In this paper, we investigate the bound-state properties of ECSC potential using a new 
perturbative formalism [24] which has been claimed to be very powerful for solving the 
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Schrodinger equation to obtain the bound-state energies as well as the wave functions in 
Yukawa or SSC potential problem [24,25] in both bound and continuum regions. This novel 
treatment is based on the decomposition of the radial Schrodinger equation into two pieces 
having an exactly solvable part with an addiitional piece leading to either a closed analytical 
solution or approximate treatment depending on the nature of the perturbed potential. 

The contents of this paper is as follows. In section 2 we breifly outline the method 
with all necessary formulae to perform the current calculations. In section 3 we apply the 
approach to the Schrodinger equation with the ECSC potential and present the results 
obtained analytically and numerically for the bound-state energy values. Finally, in section 
4 we give our concluding remarks. 



2. The Method 



For the consideration of spherically symmetric potentials, the corresponding Schrodinger 
equation, in the bound state domain, for the radial wave function reads 

h 2 <(r) 



2m (r) 
with 

V(r) 



V(r) - E n , 



(2) 



V (r) + 



h 2 £(£+!) 



2m 



+ AV(r), 



(3) 



where AV(r) is a perturbing potential and ip n (r) = Xn{r)u n (r) is the full radial wave func- 
tion, in which Xnij') is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas u n (r) is a moderating function corresponding to the perturbing potential. 
Following the prescription of Ref. 24, we may rewrite (2) as 



+ 



+ 2- 



V(r) - E n . 



(4) 



2m\Xn(r) u n {r) Xn(r)u n (r) j 

The logarithmic derivatives of the unperturbed Xn(j') and perturbed u n {r) wave functions 
are given by 
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w (r) - and AW - ( 5 ) 



which leads to 



Wl(r) - -=W n (r) 



V (r) + 



h 2 £(£+1) 



2m r 



2 



(6) 



where e n is the eigenvalue for the exactly solvable potential of interest, and 



( <ir) 2 Xn(r)<(r) \ = AW 2, s _ J^AW'ir) + 2WJr) AWJr) = AV(r) - Ae n , 
2m \u n (r) X n(r)u n (r)J nK } ^ nK ' 1 ' 1 ' 1 ' 

(7) 

in which Ae n = + + • • • is the correction term to the energy due to AV(r) and 
E n = e n + Ae n . If Eq. (7), which is the most significant piece of the present formalism, can 
be solved analytically as in (6), then the whole problem, in Eq. (2) reduces to the following 
form 

W n (r) + AW n (r) f - JL(W n (r) + AW n (r))' = V(r) - E n , (8) 

which is a well known treatment within the frame of supersymmetric quantum theory 
(SSQT) [26]. Thus, if the whole spectrum and corresponding eigenf unctions of the un- 
perturbed interaction potential are known, then one can easily calculate the required super- 
potential W n (r) for any state of interest leading to direct computation of related corrections 
to the unperturbed energy and wave function. 

For the perturbation technique, we can split the given potential in Eq.(2) into two parts. 
The main part corresponds to a shape invariant potential, Eq. (6), for which the super- 
potential is known analytically and the remaining part is treated as a perturbation, Eq. 
(7). Therefore, it is obvious that ECSC potential can be treated using this prescription. In 
this case, the zeroth-order term corresponds to the Coulomb potential while higher-order 
terms consitute the perturbation. However, the perturbation term in its present form cannot 
be solved exactly through Eq. (7). Thus, one should expand the functions related to the 
perturbation in terms of the perturbation parameter A, 
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AV(r;A) = 5>VS(r), AW n (r; A) = ]T W«(r), ^(A) = ^A^, (9) 

j=l i=l i=l 

where i denotes the perturbation order. Substitution of the above expansions into Eq. (7) 
and equating terms with the same power of A on both sides up to 0(A 3 ) gives 

2W n (r)WM(r) - = Vi(r) - (10) 

W^{r) + 2W n {r)w£\r) - *dWpr) = _ ^ 

2(W B (r)W?>(r) + w£\r)w£\r)) - = V 3 (r) - . (12) 

Hence, unlike the other perturbation theories, Eq. (7) and its expansion, Eqs. (10-12), 
give a flexibility for the easy calculations of the perturbative corrections to energy and 
wave functions for the nth state of interest through an appropriately chosen perturbed 
superpotential. 



3. Application to the ECSC Potential 

Considering the recent interest in various power-law potentials in the literature, we work 
through the article within the frame of low screening parameter. In this case, the ECSC 
potential can be expanded in power series of the screening parameter 8 as [10] 

V(r) = - exp(-5r) cos(5r) = -(-)£ V t (5r)\ (13) 

where the perturbation coefficients Vi are given by 

V 1 = -1,V 2 = 0, V 3 = 1/3, V 4 = -1/6, V 5 = 1/30, • • • . (14) 

We now apply this approximation method to the ECSC potential with the angular momen- 
tum barrier 



V(r) = - (-) e- Sr cos(5r) + 



2mr 5 



V (r) + ^ i+1 ^ h 
2mr 2 



+ AV(r), (15) 



where the first piece is the shape invariant zeroth-order which is an exactly solvable piece 
corresponding to the unperturbed Coulomb potential with V (r) = —A/r while AV(r) = 
A5 — (A5 3 /3)r 2 + (A<5 4 /6)r 3 — (A<5 5 /30)r 4 + • • • is the perturbation term. The literature is 
rich with examples of particular solutions for such power-law potentials employed in different 
fields of physics, for recent applications see Refs. [27,28] . At this stage one may wonder why 
the series expansion is truncated at a lower order. This can be understood as follows. It is 
widely appreciated that convergence is not an important or even desirable property for series 
approximations in physical problems. Specifically, a slowly convergent approximation which 
requires many terms to achieve reasonable accuracy is much less valuable than the divergent 
series which gives accurate answers in a few terms. This is clearly the case for the ECSC 
problem [29]. However, it is worthwhile to note that the main contributions come from the 
first three terms. Thereby, the present calculations are performed up to the second-order 
involving only these additional potential terms, which suprisingly provide highly accurate 
results for small screening parameter 5. 



A. Ground State Calculations (n = 0) 



In the light of Eq. (6), the zeroth-order calculations leading to exact solutions can be carried 
out readily by setting the ground-state superpotential and the unperturbed exact energy as 
h £+1 fm A , n , mA 2 



W n=0 (r) = 



+ 1 [TTT A 
f 2m r 



2 (£+l)h- 



"2/i 2 (n + ^ + l) 2 ' 



n = 0,1,2,.... 



(16) 



and from the literature, the corresponding normalized Coulomb bound-state wave function 
[30] 



Xn(r) = A£V +1 exp [-(3r\ x L 2 n i+1 [2(3r\ , 



(17) 



in which 



and L k n (x) = £ 



2mA 



(n+e+i)h 2 



£+1 



1 



is a normalized constant, (5 = 



mA 



f_1 \ m (n+k)\ 
m=0\ L ) (n~m)\(m+k)\m 



jX m is an associate Laguarre polynomial function [31]. 
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For the calculation of corrections to the zeroth-order energy and wavefunction, one needs 
to consider the expressions leading to the first- and second-order perturbation given by Eqs. 
(10-11). Multiplication of each term in these equations by Xn( r )> an d bearing in mind the 
superpotentials given in Eq. (5), one can obtain the straightforward expressions for the 
first-order correction to the energy and its superpotential: 



= r Xl(r) 

J —CO 



A5 3 



r 2 dr, (r) 



/2m 1 



xl{x) 



A5 3 



-x 



dx, (18) 



and also for the second-order correction and its superpotential: 



= Xl(r) 

J — CO 
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-r 3 - W^ 2 (r) 



dr, 



Wf ) (r) 



^2m 1 



/r 
Xlix) 



dx , 



(19) 



for any state of interest. The above expressions calculate W^\r) and W^\r) explicitly 
from the energy corrections E^ and respectively, which are in turn used to calculate 
the moderating wave function u n (r). 

Thus, through the use of Eqs. (18) and (19), after some lengthy and tedious integrals, 
we find the zeeroth order energy shift and their moderating superpotentials as 

E (i) ft 4 (I +1) 2 (I + 2) (21 + 3) ^ 
6Am 2 



p( 2) _ h 6 (£+lf(£ + 2) (21 + 3) (21 + 5) , 4 

h 10 {l + l) 6 {£ + 2) (21 + 3) (8l 2 + 371 + 43) 
72A 4 m 5 



s 6 , 



W£\r) = 



h(£+l)5 3 r( h 2 (£+!)(£ + 2) 



3V2m 



Am 



Wg\r) = 
in which 



hS 4 cr 
2v / 2m 



5 2 r 3 + ar 2 + b 



r + 



k 2 (£ + !)(£ + 2) 
Am 



* ('+1)^(2) 



^A 



) 



(20) 



h 2 (£ + l)(3£ + 7)5 2 3Am 
a = : — — , b 



c = 



Am 
h 2 (£ + lf 
9Am 



n 2 (e + iy 



ft 4 (l+l) 2 (8l 2 + 37l + 43)5 2 _ 3 (21 + 5) 
2A 2 m 2 2 (£+1) 



(21) 



Therefore, the analytical expressions for the lowest energy and full radial wave function of 
an ECSC potential are then given by 

E n=0/ = E^o, + A5 + E^ + E^ + • • • , ip n=0 ,e(r) ~ Xn=o,e( r ) u n=oA r )i 



in which 

u„ =0 ,e( r ) ~ ex P 



f 2m r 



h 



(W^ (x) + W^ 2) (x)) dx 



Hence, the explicit form of the full wave function in (22) for the ground state is 



^n=oA r ) = 



2mA 

(£ + i)h 2 



Am 



-r m exp(P(r)), 



(£+l) 2 \ h 2 (2£ + l) 

with P(r) = J2i=iPi rt is a polynomial of fifth order having the following coefficients: 

(I + 1) (2) Am _9{£ + 2) 1 .„„ 1 

A (£+l)h 2 



(22) 



(23) 



(24) 



4(£+l) 
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8 



10 



Pi — ; £ ( f ~ ^ + A 7~- /■'- 7- ^7 c2 ^ 4 ? P3 = - cc ^ 4 > = -° C( ^ = 777 c ^ 6 > 

in which d — b + ^^_^ 2<5 and other parameters are given in (21). 



(25) 



B. Excited state calculations (n > 1) 

The calculations procedures lead to a handy recursion relations in the case of ground states, 
but becomes extremely cumbersome in the description of radial excitations when nodes 
of wavefunctions are taken into account, in particular during the higher order calculations. 
Although several attempts have been made to bypass this difficulty and improve calculations 
in dealing with excited states, (cf. e.g. [32], and the references therein) within the frame of 
supersymmetric quantum mechanics. 

Using Eqs. (5) and (17), the superpotential W n (r) which is related to the excited states 
can be readily calculated through Eqs. (18) and (19). So the first-order corrections in the 
first excited state (n = 1) are 

8 



E (i) = ^(£ + 2Y(£ + 7) (21 + 3) ^ 
6Am 2 



A, „ _ M< + 2)f f ft'« + 2)«+3) _ 2f,*(e + W + 2? \ 
1 ''' ' 3^[2~^l \ Am A 2 m 2 J ' 

Consequently, the use of the approximated first two terms of w[ l \r) in the preceeding 
equation in (19) gives the energy correction in the second-order as 

(2) ^ ft 6 (I + 2) 3 (I +11) (21 + 3) (21 + 5) 4 

hl 2+W d 

ft 10 {£ + 2) 6 (£ + 3) {21 + 3) (7l 2 + 101l + 211) j6 



5 6 , (28) 



72A 4 m 5 

or 

(2) ^ fe 6 (^ + 2) 3 (f +11) (2f + 3) (2l + 5) ^ 4 
^ 24+V^ d 

_ft 10 {£ + 2) 5 (16£ 5 + 294£ 4 + 1795£ 3 + 5085£ 2 + 6878£ + 3568) 

72A 4 m 5 

if all the terms in (26) are used. Therefore, the approximated energy value of the ECSC 
potential corresponding to the first excited state is 

E n=1/ = E[ 0) + A5 + E? + E[ 2) + ■■■. (29) 

The related radial wavefunction can be expressed in an analytical form in the light of Eqs 
(18), (19) and (22), if required. The appromation used in this work would not affect con- 
siderably the sensitivity of the calculations. On the other hand, it is found analytically that 
our investigations put forward an interesting hierarchy between W^\r) terms of different 
quantum states in the first order after circumventing the nodal difficulties elegantly, 
w(1){r) ^ n(n + £+l)5* f 2 + n 2 (n + £+l)(n + £ + 2) r 2h\n + £){n + £ + l) 2 \ 



3V2m I Am A 2 m 2 



(30) 



which, for instance, for the second excited state (n = 2) leads to the first-order correction 

E W = & (£ + S) 2 (£ + 2) (2l + 23) j3 
6Am 2 



<>(r) M _ M< + 3) f f ^ + 3)« + 4) 2ftV + 2)« + 3)» 1 



3y/2m { Am A 2 m 2 , 

Thus, the use of the approximated first two terms of W^ir) in the preceeding equation(19) 
gives the energy correction in the second-order as 

.2) h 6 (£ + 2) (£ + 3) 2 (2£ + 5) {2£ 2 + 45£ + 153) _ 4 

tl — 



2AA 2 m 3 

h 10 (£ + 2){£ + 3) 5 (16£ 4 + 474£ 3 + 3879£ 2 + 121 18£ + 12873) 



<5 6 . (32) 



72A 4 m 5 

Therefore, the approximated energy eigenvalue of the ECSC potential corresponding to the 
second excited state is 

E n=2 ,e = 4 0) +A5 + + Ef ] + ■■■. (33) 

For the numerical work, some numerical values of the perturbed energies of the Is and 2s 
states, in the atomic units we take H — m — A — 1, for different values of the screening 
parameter 5 in the range < 5 < 0.10 are displayed in Tables 1 and 2, respectively. The 
results are consistent to order 5 6 with earlier results obtained by applying different methods 
in Refs. [9,22,23]. Further, we display the results for the energy eigenvalues of 2s, 2p, 3s, 
3p, and 3d states in Tables 3 and 4. Our results are then compared with accurate energy 
eigenvalues obtained by other authors. Thus, through the comparison of our results with 
those of Refs. [9,10,22,23] for large n and £— values and small screening parameter values 
yields indeed excellent results. 

On the other hand, we take A = \[2 and 5 = \[2G. Cosequently, we compute the 
binding energies {—E n ^) of the lowest-lying states, Is to 3d, for various values of S. Hence, 
the detailed analysis of the results in terms of various domains of parameters A and 5 of 
ECSC potential are displayed in Table 5. For further study of the bound-state energies and 
normalizations with analytical perturbation calculation in Table 6. We consider A = Z = 4, 
8, 16, 24 in order to cover the range of low to high atomic numbers. For low strength of 
A — Z, the energy eigenvalues nobtained are in good agreement with the other methods 
for low values of the screening parameter S. Obviously, when 8 is small the Coulomb field 
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character prevails and the method has been adjusted to that. However, the results become 
gradually worse as A and/or 5 are large. 

4. Concluding Remarks 

We have shown that the bound-state energies of the exponential cosine screened Coulomb 
(ECSC) potential for all eigenstates can be accurately determined within the framework of a 
new approximation formalism. Avoiding the disadvantages of the standard non-relativistic 
perturbation theories, the present formulae have the same simple form both for ground and 
excited states and provide, in principle, the calculation of the perturbation corrections up 
to any arbitrary order in analytical or numerical form. 

Additionally, the application of the present technique to ECSC potential is really of great 
interest leading to analytical expressions for both energy eigenvalues and wave functions. 
Comparing various energy levels with different works in the literature we find that this 
treatment is quite reliable and further analytical calculations with this non-perturbative 
scheme would be useful. In particular, the method becomes more reliable as the potential 
strength increases. 
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TABLES 

Table 1. Comparison of bound energy eigennvalues for < S < 0.1 for the Is state in atomic units. 



s 


1/N [22] 


Dynamical [9] 


Shifted 1/N [23] 


E n ,e 


0.01 


-0.490 001 


-0.490 001 




-0.490 000 9 


0.02 


-0.480 008 


-0.480 007 8 


-0.480 007 83 


-0.480 007 8 


0.03 


-0.470 026 


-0.470 026 




-0.470 025 9 


0.04 


-0.460 061 


-0.460 060 9 


-0.460 061 01 


-0.460 060 8 


0.05 


-0.450 117 


-0.450 117 4 




-0.450 117 2 


0.06 


-0.440 200 


-0.440 200 4 


-0.440 200 57 


-0.440 200 


0.07 


-0.430 313 






-0.430 313 4 


0.08 


-0.420 461 


-0.420 463 6 


-0.420 463 86 


-0.420 461 7 


0.09 


-0.410 647 






-0.410 648 8 


0.1 


-0.400 875 


-0.400 883 9 


-0.400 884 21 


-0.400 878 5 
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Table 2. Comparison of bound energy eigennvalues for < 5 < 0.1 for the 2s state in atomic units. 



s 


1/N [22] 


Dynamical [9] 


Shifted 1/N [23] 


E n ,e 


0.01 


-0.115 013 


-0.115 013 5 




-0.115 013 4 


0.02 


-0.105 103 


-0.105 103 6 


-0.105 103 61 


-0.105 103 3 


0.03 


-0.095 334 


-0.095 336 6 




-0.095 334 6 


0.04 


-0.085 755 


-0.085 769 


-0.085 769 59 


-0.085 762 1 


0.05 


-0.076 406 


-0.076 449 7 




-0.076 432 6 


0.06 


-0.067 311 


-0.067 421 7 


-0.067 426 08 


-0.067 390 


0.07 


-0.058 482 






-0.058 680 


0.08 


-0.049 915 


-0.050 392 2 


-0.050 408 25 


-0.050 357 6 


0.09 


-0.041 598 






-0.042 494 5 


0.1 


-0.033 500 


-0.034 967 7 


-0.035 004 67 


-0.035 188 
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Table 3. Energy eigenvalues as a function of screening parameter 5 for the 2s and 2p states in 
atomic units. 



State 


5 


£[10,10] [10] £[10,11] [10] Pertur.[10] 


Variational [10] Shifted [23] 




E n ,e 




9« 


n 1 n 

u. ±yj 


-0.034 941 


-0.034 941 


-0.034 425 


-0.034 935 


-0.035 004 


fi7 


— n rn^i 

U.UOtJ 


lOO u 


9n 
zp 




-0.032 469 


-0.032 469 


-0.032 042 




-0.032 470 


1 5 


— n rn9 


u i o o 


9s 


n 08 


-0.050 387 


-0.050 387 


-0.050 222 


-0.050 384 


-0.050 408 


95 


—0 050 


357 fi 


9n 
zp 




-0.048 997 


-0.048 997 






-0.048 996 




U.U4:0 


QQQ Q 


2s 


0.06 


-0.067 421 


-0.067 421 


-0.067 385 


-0.067 421 


-0.067 426 


08 


-0.067 


390 


2p 




-0.066 778 


-0.066 778 






-0.066 777 


29 


-0.066 


761 1 


2s 


0.04 


-0.085 769 


-0.085 769 


-0.085 767 


-0.085 769 


-0.085 769 


59 


-0.085 


762 1 


2p 




-0.085 591 


-0.085 591 






-0.085 559 


13 


-0.085 


552 


2s 


0.02 


-0.105 104 


-0.105 104 


-0.105 104 


-0.105 104 


-0.105 103 


61 


-0.105 


103 3 


2p 




-0.105 075 


-0.105 075 


-0.105 075 




-0.105 074 


64 


-0.105 


074 4 
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Table 4. Energy eigenvalues as a function of screening parameter 5 for the 3s, 3p and 3d states in 
atomic units. 



State 5 £[10, 10] [10] £[10, 11] [10] Pertur. [10] Variational [10] Shifted [23] 



3s 


0.06 


-0.005 461 


-0.005 462 


-0.004 538 -0.005 454 


-0.005 666 38 


-0.007 077 8 


3p 




-0.004 471 


-0.004 472 




-0.004 492 33 


-0.005 405 8 


3d 




-0.002 308 


-0.002 309 




-0.002 313 56 


-0.002 924 


3s 


0.05 


-0.011 576 


-0.011 576 




-0.011 685 44 


-0.011 952 3 


3p 




-0.010 929 


-0.010 929 


-0.010 538 


-0.010 939 85 


-0.011 111 7 


3d 




-0.009 555 


-0.009 555 


-0.009 292 


-0.009 555 42 


-0.009 694 


3s 


0.04 


-0.018 823 


-0.018 823 


-0.018 707 -0.018 822 


-0.018 867 16 


-0.018 858 6 


3p 




-0.018 453 


-0.018 453 




-0.018 457 05 


-0.018 450 5 


3d 




-0.017 682 


-0.017 682 




-0.017 682 08 


-0.017 691 


3s 


0.02 


-0.036 025 


-0.036 025 


-0.036 022 -0.036 025 


-0.036 027 38 


-0.036 021 3 


3p 




-0.035 968 


-0.035 968 


-0.035 965 


-0.035 967 71 


-0.035 964 


3d 




-0.035 851 


-0.035 851 


-0.035 849 


-0.035 850 66 


-0.035 849 


Table 5. Energy eigenvalues of the ECSC potential in units of h = 


m = 1, A = 2 1 / 2 


and 5 = GA. 


G 


State 


— £o,o 


State — £i,o 


State — £o,i State 


—Ei t \ State — £q,2 


0.002 


Is 


0.996 000 2s 0.246 000 2 2p 0.246 000 1 3p 


0.107 112 3d 


0.107 111 4 


0.005 




0.990 000 2 


0.240 003 4 0.240 002 4 


0.101 125 5 


0.101 116 


0.010 




0.980 001 9 


0.230 026 9 0.230 019 3 


0.091 221 7 


0.091 147 5 


0.020 




0.960 015 6 


0.210 206 6 0.210 148 9 


0.071 928 1 


0.071 361 7 


0.025 




0.950 030 2 


0.200 395 3 0.200 285 7 


0.062 648 5 


0.061 566 5 


0.050 




0.900 234 4 


0.152 865 2 0.152 099 1 


0.022 223 5 


0.014 137 4 
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Table 6. Energy eigenvalues of the ECSC potential for all states in units of h = 2m = 1, and 
5 = 0.2 fm~ l . 



A 


t 


n 


— E n / 


A 


t 


n 


— -En,^ 


4 








3.207 029 


16 





1 


12.825 303 


8 








14.403 752 







2 


4.023 139 




1 





2.433 587 




1 


1 


4.009 505 


16 








60.801 938 


24 





1 


31.217 455 




1 





12.818 287 







2 


11.279 786 


24 








139.201 31 




1 


1 


11.269 899 




1 





31.212 563 




1 


2 


4.412 177 




2 





11.249 961 




2 


1 


4.380 887 












2 


2 


1.411 568 
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